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Abstract 

In this paper I will approach the computation of the maximum density of regular lattices 
in large dimensions using a statistical mechanics approach. The starting point will be some 
theorems of Roger, which are virtually unknown in the community of physicists. Using his 
approach one can see that there are many similarities (and differences) with the problem of 
computing the entropy of a liquid of perfect spheres. The relation between the two problems is 
investigated in details. Some conjectures are presented, that need further investigation in order 
to check their consistency. 

1 Introduction 

The aim of this paper is to study a well known and celebrated problem, i.e. which is the 
maximum density of hard spheres when they are packed on a regular lattice. In two and 
three dimensions the solution is well known; the lattices with maximal packing density are the 
hexagonal and the ffc lattices respectively. 

In generic dimensions the result for the maximal packing density in not known [l] : an lower 
bound on the maximal packing density has been established by Minkowski [2] and it has only 
marginally improved by later studies. 

In this paper, after a brief review of the established mathematical results, we will show how 
the usual techniques of statistical mechanics may be used in this problem. We will not be able to 
find the maximal packing lattice density in high dimensions, but we believe that our work may 
be an useful step in this direction, especially for pointing the connections between this problem 
and the statistical mechanics disordered systems, like spin glasses and structural glasses. As we 
shall see later, there are very interesting relations among this problem and the thermodynamic 
properties of an hard sphere liquid. 

In section II we shall introduce the basic concepts and definitions which will be used later, in 
particular we will introduce the Roger's measure over the space of all possible lattices. In section 
III we will present the main mathematical results, due to Roger in particular we will show 
how to compute the moments of a function of the lattice using Roger's measure. In section IV 
we perform some simple computations in order to become more familiar with Roger's results. 
In section V we show how the whole approach can be simplified if we introduce connected 
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correlations functions. In section VI we present a crucial technical conjecture that is needed 
to make further progresses. Finally in the last section we start a comparison between the the 
packing problem and the hard sphere liquid; we present a possible conjecture for the relation 
between these two problems in the the infinite dimensional limit. Further work is needed to 
find out if the conjecture is consistent. Finally we present our conclusions and some comments 
on possible developments. An appendix is devoted to the computation of some integrals in the 
infinite dimensional limit. 

2 Some definitions 

Let us start by defining the problem and by establishing our notations. The first choice is how 
to parametrize the regular lattices having cells of unit volume. 

This can be done by considering the set of unimodular real square matrices Dx D (i.e. those 
matrices with determinant equal to 1). To each matrix A of this set we can associate a regular 
lattice, that is given by all the points of the form 

Xi = Yl K,knk = (An)i, (1) 
k 

where the Uk take all the possible integer values (both positive and negative). This correspon- 
dence is not one to one: the same lattice can be obtained by infinitely many different matrices. 

The condition of unimodularity (i.e. det(A) = 1) implies that the volume of the fundamental 
cell of the lattice is also equal to one. 

If we consider non overlapping spheres centered on the points of the lattices, the maximum 
allowed diameter (TZ{A)) is given by 

7^(A)^ = min|Anp = min^ nj^j^ytn^, (2) 

i,k 

where the positive matrix A is given by 

^i,fc = X! (3) 

3 

and the minimum is done over all the possible choice of the integer D-dimensional vector n, 
with the exclusion of the origin. Indeed the quantity defined in eq. ([2]) is the minimum distance 
of a point of the lattice from the origin. 
Our aim is to compute 

7^M = max7^(A). (4) 

We remark that the computation of the minimum in equation ([2j) for generic A or ^ is an 
NP-hard problem for large dimensions D. Of course there are some matrices (e.g. diagonal ones) 
for which the computation is quite simple. On the other end, in cases where the off diagonal 
elements are large, the computation may becomes quite complex. This fact may suggest that 
some of the techniques used in spin glasses (where one studies the statistical mechanics of some 
NP complete problems [3]) may be relevant also in this case. 

In order to analyze better this problem we introduce the notation 

(/(ki))A = E/(i^wi) = E/(i^^i)' (5) 

n n 

where the sum is done over origin excluded. 
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The quantity TZ{A) can computed as 



n{Af = - lim ln(F(/3)A) 

= (6) 

In other words |xp is the Hamiltonian, 

Z{P)a ^ (exp(-/?|x|2))A (7) 

is the partition function and 7^(A)^ is the ground state energy. 

The typical problem in spin glasses would be the computation of the average value (over A) 
of the free energy associated to the previous defined partition function. This is not the problem 
we face here. In this language the maximum packing density problem is a minimax problem, 
i.e. it consists in finding the matrix A (or equivalently A) such the ground state energy is as 
large as possible. This problem has been recently addressed for spin glasses in ref. [5]. 

We can formulate the same problem in an alternative way if we consider the quantity 

K{R)a = {9{R-\x\))a, (8) 

i.e. the number of points, excluded the origin, inside a sphere of radius R. It is obvious that 

K{R)a = for i^<7^(A), 

K{R)a > for R> 7^(A). (9) 

Indeed the radius of the largest sphere that is centered in one point does not intersect the other 
points is the diameter of the largest spheres, that do not overlap, if they are centered on the 
points of the lattice. 

If at given A we were able to compute (•)a, we would immediately obtain the wanted result 
after a maximisation over A of the appropriate quantity. As we have already remarked the 
computation of these expectation values for a given lattice is a nasty problem that cannot 
simply solved analytically. As usual in the case of random problems we introduce an ensemble 
of problems and we try to compute the average over this ensemble. This is usually a goal that 
we can reach analytically. 

In order to formulate the statistical problem in a more precise way we must firstly a measure 
(i/u(A) on the space of all possible lattices; we will specify later the form of this measure. We 
define 

fW) = J dfi{A)F{A). (10) 

If the measure dfi{A) does not vanish near the maximum of TZ{A), we can extract the value 
of TZm from the properties of the appropriate averages. We have two different possibilities to 
perform this last step: 

• We can use the formula 

7^i, = -lim lim (11) 

We need to compute the average partition function of n replicas of the models in the slightly 
unusual limit n — > — oo. Although there are techniques to perform this computation [5j it 
would be more natural to investigate if the r.h.s. of this last equation can be computed in 
the limit where the dimension D goes to infinity using the standard replica method that 
works very well in the case n = 0. 
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• We can alternatively compute the moments of K{R)\, this quantity being defined in eq. 

m-- 

K^''>=K{R)%= J2 P{k,R)k' . (12) 

/c=0,oo 

The function P{k, R) is the probability for a lattice (randomly chosen with the measure 
dfi{A)) of having k points inside a sphere of radius R. As far as the origin (n = 0) is 
excluded, the function P{k, R) is different from zero only for even k (if a point of the 
lattice is inside a sphere also it opposite is inside the same sphere). When R> IZm there 
are no lattices with no points inside a sphere of radius R. It is evident that 

P(0,i?)>0 for RkTZm, 

P{0,R)=0 for R>1Zm. (13) 

The program consists in reconstructing the function P{k^ R) from its moments: in this 
way one can find the value of TZm- 

Before going on with these computations, we must chose the measure dfi{A). Three possi- 
bilities come to our mind. All of them are reasonable, however, depending on the technique we 
use, for some choices the computation of the average over A may be easier. 

• There is a natural definition (due to Siegel [7]) of the measure over unimodular matrices, 
restricted to the fundamental region, the fundamental region being defined in such a 
way that each lattice may be represented in one and only one way by a matrix in the 
fundamental region. Computations with this measure are technically rather difficult and 
they will not considered here. 

• We can consider the set of matrices A (introduced by Roger) that at fixed uj depend on 
D — 1 a-variables. The action of such a matrix on a vector is defined as 

(An)j = ujTii for i < D, 
{An)D = r] inD+ atrn , (14) 

where the unimodularity condition implies that 

oj^-\ = l. (15) 

The a;-dependent measure is obtained by taking a fiat measure in the interval 0-1 for each 
for each of the D — 1 a-variables. Eventually the limit u; — > is taken. It is already 
non-trivial to prove that each lattice may be represented under the previous form. The 
proofs of this and other difficult points can be found on the paper of Roger [1] and they 
will not reproduced here. 

The computations with the Roger measure are much simpler (many results are known). 
Luckily enough the two measure are equivalent for our purposes: the quantities 

im, (16) 

that will play a crucial role in our study, are the same if evaluated with the Roger measure 
or with the Siegel measure. 

• As suggested by Kurchan and Mezard we can consider the Gaussian measure. 

d^(A) oc n dAi^k exp(-7T¥(AA*))(5(det(A) - 1). (17) 

i,k 
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where 7 is an arbitrary parameter. In this case the some of the computations can be done 
with the same techniques used in many physical problems, e.g. spin glasses [3]. The same 
measure may be useful in numerical simulations. 

As we shall see later, a simple computation shows that in the case of the Roger measure we 
have that 

Tr)=ld^xf{x). (18) 



If we apply this result to the case /(x) = 9{R — \x\) we find that 

Xp)= P{k,R)k = Vn{R) , (19) 

A;=0,oo 

where Vi:){R) = R^it^/'^T{D /2)^^ is the volume of the D-dimensional sphere of radius R, that 
obviously depends also on D. If we remember that the function S{R)/^ may take only even 
integer values, we find that P(0, R) must be different from zero in the region where 

Vd{R) < 2, (20) 

that is a consequence of the celebrated Minkowski theorem in large dimensions 0. 
The proof is rather simple. Indeed using eq. (llSp we have that: 

Vd{R) = K{R)= E P(.k,R)k>2 ^ P{k,R) = 2{l-P{0)). (21) 

k=2,oo fc=2,oo 

Let us call Rc the Z)-dependent value of R such that 

Vd{Rc) = 1, (22) 

its value being given by the condition 

RPtt^/^ 



T{D/2) 



1 . (23) 



It is convenient to measure all the quantities in units of Rc when the dimension D go to infinity. 
At this end we define: 

r = R/Rc, TM = Um/Rc . (24) 

The problem we face consist in finding the limit to infinite dimensions D of that is 
supposed to exist. The following bounds are known for large D: 

1 < < 1.322. (25) 

The lower bound is the Minkowski theorem while the upper bound (the Kabatiansky- Leven- 
shtein bound), comes from a totally different approach that we cannot discuss in details for 
reasons of space [8j. A simpler upper bound has derived by Levenshtein ^: 

TM <^ = 1-3591 . (26) 

This bound is weaker of the previous, but the proof is more direct. It is still simpler to prove 
the Roger bound ^: 

rM<V2 = 1.414 . (27) 



more precise consequence of the Minkowski theorem is that P{Q, R) is different from zero when C,{D)Vd{R) < 2 
where C(^) is the Riemann zeta function (C(oo) = 1). As a consequence we have that TZm > Rh where C(-D)Vd(-Ri3) = 
2. 
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3 Some known results 



3.1 Roger's main theorem 

The aim of this section is to recall some known results of the value of /'^ using the Roger measure 
in the case where n < D. We shall see later how this annoying constraint (n < D) may be 
removed. 

Following eq. (|14|) let us write An = X as using equation p4|) . Using equation eq. p8|) we 
find that 

/oo 
d'^xfix) . (28) 

We segue into the computations of the higher moments of /. Fortunately the appropriate 
computations have been done in rigorously way by Roger. One finds that his results are relatively 
simpler in position space. 

The strategy for controlling the previous formula consists of a few steps: 

• We classify all sets of s vectors n according to their linear dependence. 

• We perform the average and the sums inside each class. 

• At the end we write the sum over all possible classes. 
The crucial theorem is based on the following two lemmas. 

• One can prove 

n (E/(^^')) = ( / d^xfix)) ' , (29) 

where the sum is done over all sets of k vectors n'^ belonging to that are linearly 
independent (as usual the origin, i.e. n'^ = 0, never appears). The formula is valid only 
for s < D. (It is obvious that we cannot find more than D linearly independent vectors in 
dimensions D.) This formula has a very simple meaning. Each of the s points n'' may be 
carried in any point of the space independently on the other one provided that they are 
linearly independent. 

• In order to obtain the final result need also to consider the average of similar sums restricted 
to the case where the vectors n are linearly dependent. 

More precisely we introduce s vectors n that are linear combination of s — h linear inde- 
pendent vectors and their linear dependence is specified by a matrix M. In other words 
the s vectors n span a (s — h)D dimensional space and they satisfy the h different linear 
conditions: 

Mj^kn\ = for j = 1, /i , (30) 

k=l,s 

where the matrix M has integer elements and it is irreducible (there is no integer matrix 
M' such that M' = pM with p integer) . In this case we find that 



n E/(A^') =^(M) f n d^x'^fix") n ( E > (31) 

k=l,s \ „fc / k=l,s j=l,h \fc=l,m / 

where the sum is done on vectors that satisfy the condition (|3Up (that depend on the matrix 
M and have s — h linear independent components). The quantity is a normalization 
factor that is equal to 1 in many cases and that for our purposes may be taken equal to 1. 
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We now we all the elements to implement the strategy to compute the moments. We must 
consider all possible linear dependence of the vectors n and we must reduce the sum over all 
the values of n to a sum over all possible linear dependencies. 
We can now quote the main theorem of Roger [3] : 

<7^= E E-^(^) / n d'^x'^fix^) n '^'^ I E M^^kA , (32) 

h=0,s-l M k=l,s j=l,h \k=l,s J 

where the sum is done over all the sets oi h x s matrices M corresponding to different linear 
conditions El- 



3.2 The case s = 2 

Let us show how this strategy works in the case s = 2. Here we have to compute 



E . (33) 

ni,n2 

Now we have two possibilities: 

• The two vectors n are linearly independent. We obtain the following contribution 

jD^ rD 



d^Xid^X2/(xi)/(x2) . (34) 

We have to consider the case where the two vectors n are linearly dependent. In this case 
we can write the constraint in an unique way as 

qin^ + q2n^ = 0, (35) 

with if we restrict ourselves to the case of positive qi and {qi, ^2) = 1 (i-e. the pair qi and 
q2 is irreducible) H. We finally finds the following contribution 



d^x(i^y/(x)/(2/)(5^((7ix + q2y) = J d"" x f (q^x) f {q2x) . (36) 
Putting everything together we find the final expression 

E /K)/(^2) = / d^xid^X2/(xi)/(x2) + E / d''xf{qix)f{q2x), (37) 

where the sum is restricted over the irreducible pairs with qi positive. 

We have already seen how this construction works for s = 2, in the next subsection we will 
consider in details the case s = 3. 



^In order to verify if two matrices M correspond to different linear conditions, Roger writes the matrices M in a 
canonical form and presents his main theorem in this form. However we do not need his theorem under this form. 
•^The expression (a, b) denotes as usually the maximum common divisor of a and b. 
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3.3 The case s = 3 

If wc apply the previous formula in this ease we have only three possibilities: h = 0, h = 1 and 
h = 2 (h being the number of linear constraints). 

• h = In this case the matrix M does not exist and the corresponding contribution is 
simply given by 

J d^a;id^a;2d^X3/(xi)/(x2)/(x3) . (38) 

• h = 1 In this case we have a linear constraint of the form 

qixi + q2X2 + qaxs = . (39) 

We have two possibilities. 

— One of the three q is equal to zero. Here we have three equal contribution. If we set 
qs = 0, we must impose that qi and q2 have no common factors (i.e. (gi, 92) = !)• We 
finally find a contribution equal to 

3 E / d''xifiqixi)fiq2x,) [ d'^xsfixs) . (40) 

— All the q are non zero. The previous condition requires that qi, q2 and qs have no 
common factors. We finally find 

E ,3-^ / d^x^f{x,) [ d^X2f{x2)f (^ia±^) . (41) 

• /i = 2 In this case we can write the two linear constraints as 

qix\ = q2X2 and gixi = 5-2X3- (42) 
We finally finds that the contribution from this case is given by 



E (9252)-^ / d^xf{x)f (^) / 



9l>92,9l,92 



9lX 
92 



(43) 



The final result is given by the sum of all the previous contributions. 

In the same way we can write explicit formulae that contain more and more terms when s 
increases. 



4 Some simple computations 

In this section we will present some simple computations, in order to familiarizes ourselves with 
the previous results. 

If we look to the previous formulae it is evident that there is an obvious prefactor that goes 
to zero very fast when D ^ 00 for many choices of the function /. Let us see what happens in 
a some simple examples and/or doing some simple approximations. 
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4.1 The Gaussian case 

Let us consider the case where the function f{x) is given by 

/(x) = (27r)-^/2exp(-^xV2). (44) 
Here the integrals can be easily done and we finds that 

< f >= r ^ + E ill + gi)-^/' (45) 

91.92 
91.92 

E {ql + il + iir^+p-""' E {ihl+qhl+qhl)-'"' 

91.92.93 91.92.91.92 



For large D we find: 



D 

f2 



21n(</>) 



^M4/^«.ma.(/3-,(2/3)-) 



D 

P3 



(46) 



^^^^ max(/3-^ (2/3^)-\ (3/?)-^) 



The problems connected with the limit D ^ oo are clear from the previous formulae. De- 
pending from the vahic of /3 different terms are the leading ones (the case /3 > 1 and (5 <2 are 
quite different). Sometimes it happens that the leading term from one contribution is smaller 
of the subleading terms of other contribution. The terms with the q's and the g^s equal to 1 are 
always the leading ones inside a given class. 

4.2 The theta function 

We consider here the case of the function: 

f{x) = e{R-x). (47) 

Our plan is to compute the leading term of each moment of / when the dimension goes 
to infinity, resuming the resulting series and extracting information on the probability P{k,d). 
The relation /^(x) = f{x) will be useful in simplifying the result. 

Let us use the previous formulae for computing the first three moments. We find in the 
same way as before that 

I 

= = VD{Rf + Vd{R) E (max(gi, \q2\r'' (48) 

91.92 

i^(3) = 77^3 = 

Vn{Rf + Vn{R)j:^^{qMr^ + Vu{R? E (^iil^^Z^ll^MM)-^ + 

91,92 qi, 92, 91, 92 ^1^2 

E (93)"-° / d^'xi f (fx2{e{zi - R)e{x2 - R)9i\qixi + 52X2! - qsR) 
91,92,93 
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We have terms that scale as different powers of V£){R). Inside each class the leading terms 
are those the variable q and g are equal to ±1. If we keep only these terms we find the simple 
result 

=< / >= Vd{R) 
^(2) =< /2 ^= y^(^)2 + 2Vd{R) (49) 
^(3) j3 >^ 

VD{Rf + 2VDiRf + ^VoiR) d^xi I d^X2e{xi - R)9{x2 - R)9i\xi + X2\ - R) 

The last integral is exponentially small for large D with respect to the other terms (see also 
the appendix): it is the probability of havingthree overlapping spheres of the same radius. This 
crucial results follow from two general facts □ that are valid when D ^ oo: 

• The measure on a D dimensional sphere is concentrated on its surface. 

• Two generic vectors (of length 1) have vanishing scalar product; therefore the sum of two 
generic vectors of length less or equal to R is at most 2^^^R. 

Also in this case (depending if R is larger or smaller that Rc), different terms dominates the 
result for each moment when the dimension becomes large. 



4.3 Summing the leading terms 

Our aim it to obtain a closed formula for the leading terms for each moment of the function 
/ in order to extract the asymptotic behaviour. Unfortunately controlling the leading term in 
each moment does not imply that we control the leading term in the probability distribution. 

From the analysis of the first three moments we conclude that in general the terms that are 
leading when the dimension goes to infinity are those that come from h linear conditions of the 
form 

Xi = ±Xj. (50) 

If we select these contributions we find that the result can be given in terms of h integrals 
of the same function. In the case of a generic even function (i.e. f{—x) = f{x)) we get: 



h=l,siyi=l,s 1=1, s 



dx^f{xY' (51) 



where the sum is done over all the sets of h integers z/j such that 

i 

and C{s,h,v) is a crucial combinatorial factor that is equal h\^^ times the number of ways in 
which we can divide s objects in h group of Vi elements (we have to specify which variables are 
equal). This number is given by 

C{s,h,v) = -—^ (53) 

For example in the case s = 3 and h = 2, we two equal contributions 

Ui = 2 U2 = 1 

V2 = 2 VI = I (54) 
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The argument is well known and plays a crucial role in the evaluation of the leading term of the Mayer expansion 
for an hard sphere gas in infinite dimensions (see [H])- 



10 



that correspond to the same term (the number of contributions compensate the term m! ^ = 
1/2). The remaining factor 3 correspond to the three contributions: 

dxidx2dx3[6{x2 - xi) + 6{xs - xi) + 6{xs - X2)] ■ (55) 

The formulae simplify if we compute the generating function of the moments. We find after 
some simple algebra that 

expy < / > PS exp(^ J dx(exp(2y/(x)) - 1)) . (56) 

If we apply this formula to the case of the function in eq. (j47p we find that 



expyK{R)A « exp (^-M^exp(2y - 1)) 



exp 



Vd{R)\ V- 1 fVDiR)\ 



k 



We finally find the simple result: 

,„,.,..,(zm)>(m)', 

while P{2k + 1,-R) =0 (The number of points of a lattice inside a sphere is alway even!). 
As a byproduct of this Poisson distribution we have that 

P(0,i?)^exp(^-^) . (59) 

Let us find the consequences of this approximate result. We have called Rc the radius of 
a sphere of unit volume. If we take a large value of D at fixed ratio r = {R/Rc) we find that 
V{rR(^ goes to zero or to infinity depending if r is smaller or greater than 1. 

• If r < 1 we have that P(0, rRc) goes to 1 apart from exponentially small corrections. With 
probability going to 1 lattices do not have points at distance smaller that Rc- 

• In the interesting case r > 1 this computation give an exponentially small but non van- 
ishing contribution for P{0,rRc). 

If we stick to this result we obtain that P(0, rRc) is always non zero and that we can find always 
a lattice that contains no points at distance less than rRc- This result cannot be correct because 
it can be proved that no such lattice exist for large D as soon r > 1.32. 

This failure is due to the fact that we have neglected exponentially small terms. This is 
allowed only if these terms do not compensate the exponential small term we have obtained for 
the probability distribution and produce a net zero result. In other words we cannot anymore 
neglect the subleading terms if the sum of the of the leading terms becomes smaller of each 
individual subleading term. It was proved by Roger that this does not happens in the relatively 
smaller region where < D {a similar region appears in the analysis of [ISj). 

Up to now we have transcribed the results of Roger in a slightly different language skipping 
the proof of his theorems. Our aim is to extend our command of the expression for P{k, R) by 
considering more terms in such a way to be able to control better the remaining integrals. This 
will be done in the next sections. 
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4.4 The appearance of the random energy model 

Before going to more precise computations, we it is convenient to study of also the partition 
function if we retain only the the terms considered in the previous subsection. We recall the 
definition of the partition function given in eq. ([7]). 

Z(/5)a = (exp(-/?|x|2))A . (60) 

If we apply the formulae of the previous subsections we find after some algebra that 



Zif3y = 2^ E E Cis,h,n) n 2-M|^)"''^' , (61) 

h=l,sni=l,s i=l,h ^'"^ * 



where (as usual) the sum is done over all the sets of h integers rii such that J2i n-i = s. 

The previous formula does not look particularly illuminating, however a close look shows 
that this formula is well known to people studying the random energy model (REM) [lO] . In 
this way we can rewrite it in a more illuminating way. At this end let us consider a new model 
where the partition function can be written as 

Zr{P) = E ewi-PEk) , (62) 

k=l,N 

where the Ek are random independent quantities distributed with the probability distribution 
P^^\E). Standard arguments shows that 

^= E E C{s,m,n) n pkV, (63) 

m=l,sni=l,s i=l,m 

where 



p\l^> = '—^ y->^-^y ■"^-> = y dEp^J^i(^) exp(-n/3^) (64) 

and we have defined 

V'^\E) = ?^. (65) 

In other words p^^\E)dE is the probability of finding an energy in the interval [E,E + dE]. 

The limit going to infinity can be done without difficulties if the function |?(^) (i?) depends 
on A^ in such a way that this function remains finite in this limit, i .e. if the following limit exists. 

p{E) = lim p^^\E) (66) 

N—*oo 

The normalization condition (J p^^\E)dE = N) implies that 

^ P{E) = oc, (67) 
In our case we have to take 

p{E) = lim p(^)(i?) = (68) 

It is immediate that within the approximation of summing the leading terms the partition 
of our model coincide (if we neglect the factor 2** in eq. (jOip □ ) with a random energy model 
where 

Pn = 2-1 i^—j =-j d^a;exp(-n/3xV2) (69) 



^The factor 2^ , that is irrelevant for most conclusions may be easily obtained by adding a fixed degeneracy ( 
of the levels. 
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and therefore 

p{E) = ^J d^x6iE - . (70) 

In other words we consider points in the D-dimensional space that are uniformly 
randomly distributed with density 1/2 inside a sphere of an appropriate radius going to infinity 
with N. The partition function is just the limit ^ oo of 

ZM = J2 exp(-/3xf,/2) . (71) 

k=l,N 

This implies that the probability distribution of the quantities |Anp correspond to a flat un- 
correlated distribution of the points An. The probability of finding an arbitrary large sphere of 
the space with no points, is always non zero although it is very small for a quite large sphere, 
this is perfectly consistent with the results of the previous subsection, where it was shown that 
there is no upper limit for the radius. 



4.5 A detailed analysis of the Random Energy Model 

With this identification we can proceed to a direct analysis of the REM. Within this approxima- 
tion we would like to follow the strategy outlined in section II, where we connect the maximal 
radius with the behaviour of Z{(3Y in the double limit /? — > oo and s —>■ — oo. 

The first method we could use consists in computing the moments of the partition function 
in the REM replica method. In the case of the REM the replica method is well established in 
this case and it gives the same results of more conventional methods. 

Let us firstly follow this route. In evaluating the l.h.s. of eq. (|6ip we consider the contribu- 
tion where all the rii are equal to u. In this case their number is k = s/u. This contribution 
gives 



Z{py « 2*s!/(M!)'=exp ( -s 



L>/21n(/3n/(27r) +ln(2) 



u 



= rsl/iul)'^ l^/rf xeM-uPx /2)^ ^^^^ 

The previous formula can be used as input to the replica method. After some computations 
(that we skip) we can try to use the appropriate saddle point method in the limit D oo. 
Within this approximation, neglecting those prefactors that have a constant limit when D ^ oo, 
we find that in the low temperature region (i.e. high beta) the moment of the partition function 
are given by: 

Yipy ^ exp{-ADsp), (73) 

where the constant A is given by 

A = Tre (74) 

In other words we find that all the realization give the same partition function with proba- 
bility one. Unfortunately we are interested to the behavior the partition function in the tail of 
the probability distribution and this require the study of the limit of large /3 and —s. In this 
situation the prefactors cannot be neglected and the analysis becomes more complicated. 

If we look to the formulae of the REM we find that for sufficiently high (3 we have that the 
partition function is dominated by the ground state and that the probability distribution of the 
minimal energy is given by by 



V{E) = p{E) exp(- / de P{e)) oc E^/^ exp{-VD/2E^/'^) 

Je<E 



(75) 
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We thus find that 



(Z/2)^ w J dEV{E)exp{-l3E) oc 

J d^xexp{-(3sx'^/2 - ^Vdx^) oc J dR^ expi-^VnR^ - ^s(3R^) (76) 

If D goes to infinity first we find the previous result. On the other hand, if we have a large 
value of (3 for negative s at fixed dimension, neglecting the prefactors we find that 

ZpF«exp(-i?(Z))(/3s)i5^) . (77) 

where 

B{D) = ^ {DVd) ^ ^D-.oo 7^e (78) 

The problems arising in exchanging the two limits D and (3 to infinity for negative s are 
clear from the previous formulae. Moreover within our approximation the limit 



Um (79) 

s— >— oo s 

does not exist, contrary to our expectations. This negative result reflect our lack of command 
of the tail of the partition function within our approximation. 

These computations clearly shows that the approximation of keeping only the leading term 
is not enough to find the most compact regular lattice and that a more refined approximation is 
needed. This will be the subject of the next sections, where we are going to present new results 
and conjectures (up to now we have only translated in a different language the results of Roger 



5 The importance of being connected 
5.1 The introduction of connected moments 

Before considering further contributions we need simplify the previous computation. At this 
end it is important to introduce the connected moments and to do some general consideration. 
Let us consider an integer valued function / and its moments defined as 

<f>=Y,P{k)k^^fs, (80) 

k 

where P{k) is the probability that the function / is equal to k. 
It is usual to introduce the generating function 

G(^)= E /«f = E E P{k)^-^= E P{k)eMzk)=eMC{z)), (81) 

s=0,oo ■ s=0,oo fc=0,oo ' fc=0,co 

where C{z) is the generating function of the connected moments: 

C{z)= J2 <r>cz'- (82) 

s=l,oo 

Let us assume that we can write 

C{z) = B ^ p{k){eMzk)-l), (83) 

fe=l,oo 
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where the quantities p{k) are non-negative and such that 

E Pi'') = 1- (84) 

A:=0,oo 

Using the theorem of composed probabihty we readily find that 

P{k) = eM-B) E ^ n E Pi''^) p( E - ^) (85) 

n=0,oo ' 1=1, n \/ci=l,oo / i=l,n 

In other words the distribution P{k) can be obtained in the following procedure: 

• We first extract a number n with a Poisson distribution with average B. 

• We extract n independent numbers ki with probability = p{k). 

• The quantity k is given by J2i=i n ^j- 
We finally find that 

P(0) = exp(-^) = exp(- lim C{z)). (86) 

We see that the problem of computing the value of P{0) is reduced to the evaluation of the 
generating function of the connected moments in a particular limit. 

On the other hand, let us assume that there exist an analytic function /c(s) such that for 
positive integer s 

f,{s) =< f >e, (87) 

and the function fc{s) does not have a nasty behavior at infinity in the complex plane. If this 
function does not have singularities in the region 5Re s > one can prove that 

lim C7(z) = /e(0). (88) 

2— » — oo 

Under the previous assumptions the task of computing P{0) is reduced to the problem of finding 
an analytic expression for the function /c(s) and of computing it at s = 0. Our strategy will 
consists in summing subsets of all the possible contributions and to evaluate the result at s = 

E. 

If we use the computation of the moments done in the previous section, where only the lead- 
ing terms where taken into account, we find the very simple result for the connected moments: 

fc{s) = J d'^xfixY . (89) 

If we take the function / equal to K{x) = 6{R — x) we find that 

/(s) = J d'^xKixY = Vd (90) 

that is obviously an analytic function of s in the whole complex plane. 



^The reader may be puzzled by this result that has a strong replica flavour. 
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5.2 How to compute the connected moments 



It is evident from the previous discussion that the direct computation of the connected moments 
is extremely interesting. The final expression for the connected moments certainly contains less 
terms than the expression for the moments: this is not a surprise to people working in statistical 
mechanics. 

We discuss now how this task may be achieved. Let us consider what happens up to the 
third moment. If we start from the formulae of the section III the same notation, we find: 

< f >c= E / d''xid''x2f{xi)f{x2)S{xiqi + X2q2) (91) 
<f>c= E d^^ld^X2d^^3f{xi)f{x2)f{x3)6{qiXi+q2X2 + q3X3)+ (92) 
^ d^xid^X2d^X3f{xi)f{x2)f{x3)6{giXi + g2X2)6{qiXi + q2X3) (93) 

91,92,31,52 

We ready see that the only terms that remains in this limit are connected (in the same way 
of connected diagrams) in the sense that they cannot be split into the products of terms that 
contain integrals that can be done independently. 

In order to see how derive a similar result for the higher moments and to stream down the 
discussion let us set all factor J\f/q to 1. In this case the general formula of section III can be 
written as 

<f>=T.f n {d''xJ{x,))G{{x}) , (94) 

G i=l,n 

where the sum is done over the appropriate set of function G{{x}). These functions are an 
appropriate product of delta functions that enforces the linear dependence of some of the x's. 

We state that a function G is connected if it goes to zero when one or more of the x go to 
infinity together. In other words we must have that 

lim G^(xi + Xci,X2 + Ac2, + Ac„) = (95) 

A— >oo 

for any set of vectors c such that at least one c is different from zero and no pairs of c are equal. 

If G is not connected it can be written as the product of its connected pieces, if now we 
take care of all the multiplicity factors and we use the standard manipulations of statistical 
mechanics we find that 

<r>c=T.J n {d''x,f{xi))G,{{x}) (96) 

Gc i=l,n 

where the sum is done only over those G that are connected (we call them Gc{{x})) 0. 



6 A crucial conjecture 
6.1 A serious difficulty 

Before going further we must face a problem that we have postponed up to now. 

^Some people may consider the previous formula nearly self evident. Indeed let us consider a function / that is the 
characteristic function of a domain of volume V. General considerations may be used to argue that < >c should 
be proportional to V and this is possible only if only connected G are present in the r.h.s of the previous equation. 
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The program we have put forward may seem witless. We have said that we can write 
equations of the form 

7^ = F{s,D), (97) 

that are vahd only for s < D (as we have seen the r.h.s of the previous equations contains 
some sums over integers and some integrals). Roger theorems (that are at the basis of our 
computations) does not give useful information for m > D for the very good reason that in this 
case both sides of equation eq. (j32l) are infinite. Indeed there are very asymmetric lattices for 
which the function < / >a is very large (e.g. there are many points very near to the origin). 
Conversely the the sum over the M is not always valid. Indeed it is easy to check that the 
formula eq. (j37p for is divergent in two dimensions as soon as /(O) 7^ 0. 

How can we write a meaningfully expression in finite dimension for exp(y < / > if we can 
compute only a few moments of the < / >? We would like to put forward a reasonable conjecture 
that should allow us to obtain an explicit form it. Before doing it, let us look better to the 
previous formulae and see what happens if we refine the previous computation by including 
some additional terms. 



6.2 A more refined computation 

Let us consider for simplicity only the connected moments. The non-connected one can be 
obtained from them by combinatorial factors. We want now to analyze in details the origin of 
the divergences that we have seen in the previous subsection. 

Let us start by considering the second moment in the case where the function is equal to 
exp(— x^). In this case we have explicit formulae. Doing the appropriate integrals on obtain 
that 

I / „ \ D/2 

As we will see in dimensions larger than two the sums are convergent. 
In the same we we obtain 

<f^>= V ( 

In dimensions three only the first term is divergent. The other term is divergent in dimensions 
two. 

This is a general feature. The most divergent term is always the one that correspond to a 
single linear constraint. We can therefore start by considering only this contribution. We want 
to arrive to a final simple expression, that can be reached in a few steps. 

Within the approximation of taking the linear constraint we have 

<7^c « E / n dx^f{^k)5{ E l^^i) (100) 

qi-Qs k=l,s i=l,s 

where by the superscript / we indicate the fact that there is no common factor among all the 
q. In our case we find 

^T^c - E (^^^)^/' (101) 

qi --qs Z^i=l,s Qi 
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We note that if H is an homogeneous function of degree —u 

Y: H{q) = ({u) Y: H{q) . (102) 
qi---qe qi—Qs 

We finally find that (within the approximation of keeping only the linear constraint) 
where 

G(i)=7r^/2 exp(-iA;2). (104) 

fe=l,oo 

Few remarks are in order. 

• In this case the function fc{s) can be written under the form 

fc{s) = I dzp{z)z% (105) 

where 

p{z) oc J jt^/^6{G{t) - z). (106) 

It is crucial to note that the expression for p(z) is well defined also in the case where some 
of the moments of the function / are divergent. The divergence of the moments of the 
function / is related to the asymptotic behaviour of the function p{z) at large z. 

• We have that for small G 



G(i)=(^-j +G{t) (107) 

where G{t) is a C°° function at i = 0"*" (it is a well known fact that can be easily proved 
using the Poisson formula). Therefore the function /c(s) has simple poles on the real axis 
at integer values of s starting from s = 3. The first pole correspond to a decrease of the 

function p(z) as z^^ at large z. 

For dimensions D > s the quantity /c(s) is an analytic function of D that has some poles 
for some values of s. Moreover /c(s) is an analytic function of both variables s and D 



6.3 A preliminary conjecture 

We can hope that these analyticity properties in D and s are true in general. If we accept 
this conjecture each given moment can be computed at sufficient high dimension, where it is 
convergent, and we can evaluate it as analytic continuation in D and ,s at the point we need it. 
The analytic continuation of the moments in the dimensions at fixed s and analytic continuation 
of the moments in the value of s at fixed dimension should coincide. If in a given dimension 
they are both singular, the form of the singularity in s and in D are related. 

If this conjecture is true, Roger formulae bring information on the function fds) also in 
the region where the definition of the moments is divergent and therefore our program is not 
witless. 

This conjecture is quite strong. It is possible to verify by explicit computations that the fds) 
function has the needed analytic properties in the region where Roger formulae are convergent. 
On the other hand it is not obvious that the function fc{s) at given dimension coincides with 
its analytic continuation from high dimensions. We do not know how such a conjecture could 
be proved in a rigorous way. It would be interesting to check if this conjecture is supported by 
a numerical study of the first moments in low dimensions. We will not attempt to further study 
this point here. 
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7 A Comparison with the gas of Hard spheres 



7.1 The hard spheres gas 

Here we try to get some tentative conclusions for the asymptotic behaviour of rM in large 
dimensions. In order to do this we have first to recall some properties of the gas of hard 
spheres. 

An object that it is quite familiar to physicist is the partition function of a gas of hard sphere 
of diameter r. 




(108) 



The partition function depends on the total number of particles N and on the volume V of the 
box where the particle are confined. Without loss of generality we can assume that V = N, i.e. 
the density is one. In the infinite volume limit we have that 

Z{V, N) ^ exp(-y 5(r)) , (109) 

where S{r) is the entropy density. The entropy is a function of the reduce diameter r and it 
should diverge when r arrives to the point (r^i) that correspond to the highest density packing 
(for r > Z = 0). It is hard problem to decide if in the infinite volume limit the highest 
density packing is a lattice packing. 

In dimensions 3 the most compact packing is a lattice packing. Therefore the partition 
function vanishes as soon as r > rM and 

lim S{r) = -oo (110) 

r-*{rM) 

A crucial phenomenon for a gas of hard spheres is the phenomenon of crystallization: there is a 
reduced radius rc such that for r > rc the configurations of the hard particles are very similar 
to those of a regular crystal. The transition form the liquid to the crystal phase is a first order 
transition. 

Generally speaking in higher dimensions we can have two possibilities 

1. The partition function of hard spheres is different from zero also in the region r > vm 
and diverges at a value r = > r^j- In this case ta is the maximal radius of a non- 
crystalline (i.e. amorphous) packing and is bigger of the maximal radius of a crystalline 
packing. This situation is quite different from the three dimensional situation. 

2. As in the three dimensional case the partition function of hard spheres becomes zero at 
rM- We can however distinguish among two cases. 

(a) There is a radius rc where the entropy has a first order transition where the system 
crystallize. 

(b) There is no crystallization transition. 

The situation become more complex if we consider the possibility of having a glass transition 
in the non-crystalline phase p3] at rc- This glass transition may be or in the metastable phase 
above rc, has happens in three dimensions, or in the stable liquid phase below rc- This glass 
phase transition is important because the virial expansion does not give information on the 
behaviour above rc- The arguments of [13] predict that rc = I- 
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7.2 A first look to the virial expansion 

The effective density of tlie particles is proportional to and excluded volume effects becomes 
larger when r increases. The usual virial expansion may be used to compute S{r) as function 
of r and of the fugacity in the gran-canonical ensemble: 

^(^' ^) = E ^(^' ^) = exp(-yF(z)) . (Ill) 

Af 

Eventually we can adjust the fugacity z in such a way to have density one. These com- 
putations are well described in the literature. Diagrams are classified according to increasing 
complexity. 

Our aim is to compare the virial expansion for the partition function with the formulae that 
we have obtained for ln(P(0, r)). Let us consider firstly the results that we obtain if we consider 
a particular class of diagrams, i.e. those of a chain. In this case one finds that 

- Sir) = V{R) + E / n Mx^^^i E ' (112) 

fc=3,oo i=l,k i=l,k 

where 

du{x) = d^xd{R - x) . (113) 

We would like to compute the quantity Siir) == ln(P(0, r)). We can follow the same arguments 
as in the section ()4.3p . One finds, among many other terms the following terns (neglecting factors 
2) 

- S{r) = V{R) + E / n e) ^( E ^^90 , (114) 

fc=2,oo i=l,k \ Qi / 1=1, fe 

where the sum over the q is done with the condition: 

• All q are different from zero. 

• Conventionally qi is positive. 

• The set of k variables qi for i = l,k is irreducible, i.e. they do not have common factors. 

• For k = 2 the terms qi = 1 = ib(/2 are absent. 

At each given k the leading term in the sum comes from those term where qi = 1 and 
= ±1. If we retain only these contributions we find (apart from factors two) exactly the same 
result that for the pressure. Therefore we find that within the approximation of considering 
only the terms coming from one linear constraint and the for each moment taking only the 
leading term, the expression for Siir) coincide with that for S{r) obtained summing a subset 
of diagrams of the virial expansion. 

This result should not be surprising and it can be simply understood: we have seen that 
the distribution of the points of a random lattice are random if they do not satisfies linear 
constraints (see section (031)) |. This implies that as soon as we consider only the probability 
distribution of a finite (not diverging with D) number of points that would be just the same of 
points that are not constrained to be on the lattice. 

®This is clearly related that the matrix A contains a very large number of parameters when tlie dimensions D goes 
to infinity and that a finite number of points may be points of the lattice 
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7.3 A naive conjectures 



Let us consider the generic contribution to the the virial expansion. Every diagrams of the virial 
expansion of S{r) with k Hues and L loops can be written be written in the form: 



where the quantities M take the value or 1. The form of the matrix M fixes the topology of 
the diagram. 

It is possible to check that all the contributions of the previous form also appears in the 
computation of SL{r) = — ln(P(0, r)), and the terms with many loops correspond to terms 
with more than one linear constraint. All the diagrams of the virial expansion appear in the 
computation of SL{r), however this last expression contains much more terms that are not 
present in the virial expansion. 

It is possible that the leading contributions in both problems are the same (as it happens 
in the case of the chain approximation) so one could may be tempted to conjecture that in the 
limit of infinite dimensions the two problems coincide and in the infinite dimensional limit we 
have that 



The previous equality would implies that in infinite dimensions the maximal density of lattice 
packing is equal to maximal density of packing without any lattice constraint. Although the 
conjecture may be true order by order in the virial expansion, we should be very careful because 
the sum of subdominant term may become dominant. I believe that this conjecture has only 
small chance to be true 

7.4 A more refined conjecture and open problems 

In order to further understand the problem we need to control the hard sphere gas in large 
dimensions D and eventually to treat the difference among the two problems as a perturbation. 
The task of computing the properties of the hard sphere gas in the infinite dimensional limit is 
non trivial (when r > 1) [12\ [T3\ I14j . In some sense we have transformed an unsolved problem 
into an other unsolved problem. The situation is not so bad, as far as the computation of the 
thermodynamical properties of an hard sphere gas may be not out of reach. One should firstly 
obtain some results in the liquid phase for r > 1, and study the properties of the glass transition 
is present. 

We recall that some progresses have been recently done by resuming the contributions that 
correspond the hypernetted chain approximation in the infinite dimensional case for r < e/2 
'|12j. Also the glass transition has been tentatively located around r = 1 [13] . 

The situation is quite complex and it may be convenient to consider three functions for large 
dimensions D: 

• S{r)) is the entropy for a gas of hard spheres as function of the radius. It is a function that 
may have singularities a various values of r. In three dimensions it has a singularity at Tc 
that corresponds to crystallization. It seem that in three dimensions at rg > Vc there is a 
glass transition, characterized by a faster increase of the pressure It is possible than in 
higher dimensions the relative position of the glass transition and the crystal transition do 

^As far as this phenomenon happens in the metastable phase, things are not so well defined; however one 
introduce slightly modified Hamiltonians that forbids crystallization [15j . 




(115) 




(116) 
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change. It is also possible that in high dimensions the crystal like transition disappears, 
but this would be quite unlikely. The value of Tc is estimated to be equal to 1 jl3J . 

• SL{r) = ln(P{0,r)D, A. e. the logarithm of the probability of finding a lattice packing 
with reduced radius r. It should diverge toward infinity when we reach the maximum 
density for a lattice packing, i.e. rjv/. The reader should notice that this function can 
be written as an integral over a finite dimensional space so for fixed D it should be a 
continuous function of the radius. Only in the infinite dimensions it could develop a real 
first order phase transition, however it could be possible to observe the premonitory signs 
of this transition by evaluating it numerically. It is also not clear if the glass transition at 
r = 1 is present in this contest. 

• The two previous quantities are not easy to compute, although they are well defined. 
One could use here a third quantity that is maybe slightly less defined, but it would be 
much easier to compute, i.e. SE{r) that is the approximate value of the entropy computed 
using integrals equations. In other words SEir) coincides with the results of the HNC 
approximation as soon the bridge diagrams are neglected. An explicit solution of the HNC 
equation has been found in in the region where r < e/2, i.e. where the Levenshtein 
bound is satisfied. 

If at some value of r the bridge diagrams become important, one could consider more 
complex equation. If only a finite number of bridge diagrams become relevant one after 
the other, this program could be implemented because if we have seen there diagrams can 
be computed analytically in very large dimensions. 

In this contest there are a few points to be investigated. 

1. The computation of SEir) is the more urgent task if we would like to get analytic predic- 
tions. While the situations is more or less clear in the HNC, it is not clear if there is a 
value of the density where bridge diagrams must be considered. 

2. One should use the analytic tools that are used to predict a glass transition in the hard 
sphere liquid, to see if they predicts also a phase transition in ^^(r). 

3. One should compute the S'£;(r). in dimensions where the computation is feasible to see if 
there are sign of a first order transition. 

If we succeed in task 1, if the questions in two have a negative answer, one can put forward the 
conjecture that 



An important test of the conjecture would be to verify is the terms that we have neglected do 
not pile up. 

In conclusion this paper has opened more problem that those it has solved, however I am 
confident that strong progress can be done on this line. The difference and the commonalities 
of the two problems we are investigating could be further clarified. 
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Appendix: on the large dimensional limit of some in- 
tegrals 

In this approach it is crucial to evaluate some integrals in the limits where the dimension D of 
the space goes to infinity. Most of the integrals we consider go to zero exponentially. This result 

is related to the fact that in large dimensions most vectors are orthogonal. Here will mainly 
discussing the exponential factor, and we will neglect the prefactors. 

Before considering the general case let us start by studying a simple example: 

A = j dn{x)dn{y)5{{x + yf - 1), (118) 

where 

d^l{x) = S{D)-^25(x^ - \)d^x (119) 

The factor 

is the surface of the D dimensional sphere and it is such that 

d^Ji{x) = 1. (121) 



I' 



If we go to polar coordinates, we introduce the angle (j) betweenx and y we find, neglecting 
prefactor that are at most powers of D 

f dcos{4>) sm{(Pf6{cos{ct)) - ^) « ( j)^/^ (122) 

The same result can be obtained if we consider the same integral with the measure 

du{x) = S{D)-^d^x9{l - x'^) (123) 

still satisfying the condition 

jdv{x) = l. (124) 

We can also consider the integral 

A = j dv{x)dv{y)S{D)-^e{{x + yf - 1). (125) 



This integral is proportional to the probability of having three sphere with each sphere touching 
each of the remaining two. 

In large dimensions the measure of a sphere is concentrated on its surface and the difference 
among a ^ or a 5 function is irrelevant. Therefore the integral have a similar limit when the 
dimension goes to infinity. 

The same result could be obtained by using the integral representation for the function 

e{z) = I — expiiaz), (126) 
J a 

where the integral is done from — oo to oo, giving a small positive imaginary part to a. In this 
way we obtain that 

A ^ V{D/2f j daida2da3ex-p{iE){iV)-^/^ = 
T{D/2f j daida2da3 exp(iS - D/2 In(iP)) , (127) 
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where 

S = ai + 02 + as, ^ = aia2 + 0203 + "sOi • 
The saddle point equations are: 

2V 



i = ^5%ti5), (128) 



plus those obtained by permutations. 

The solution to the saddle equations is given by 



D , , 

ai = 02 = as = (129) 



Substituting back we find the result eq. (jl22p . 

Alternatively, we can formally introduce the variables 

Zi = ^ for z = l,3 (130) 

that satisfy the constraint 

E = 1- (131) 

j=l,3 

We finally find 

A « r(L>/2)2 y" dY^TT^-^ eyip{iT.)dzidz2dz36{ J2 - 1)^"^''^ 

- f dy 

\D) 



i=l,3 



^Y(W~ I ^^i^^2d23'^( E - l)^"""/' ~ I dzidz2dz^5{ Y^zi- l)V-''l\ (132) 

^ ' i=l,3 1=1,3 

where 

V = Z\Z2 + -Z2-Z3 + ZzZ\ (133) 
We finally find that the saddle point is given by 

^^ = I, (134) 

and we recover the previous result. 

The same technique can be used to estimate the integrals 

An= f n du{^^)0{l - ( E ^^)') (135) 

i=l,N-l i=l,N-l 

We find 

Other integral may be estimated in the same way. We may be interested to consider the 
cluster integral: 

C = f du{xi)du{x2)du{xz)e{l - (Xi + X2f)e{l - (X2 + X^f)e{l - (X3 + Xif). (137) 



This integral is proportional to the probability of having four sphere with each sphere touching 
each of the remaining three. We find that 



C ^ r(I)/2)^ 



r(3D/2) 



/ n dz,5{ Y^z,- l)p-^/2 (138) 



i=l,6 i=l, 
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where 

V = Z1Z2Z3 + Z1Z3Z4 + Z2Z3Z4 + Z1Z2Z5 + Z1Z3Z5 + Z1Z4Z5 + ^2^4-^5 + Z3Z4Z5 + 

Z1Z2ZQ + Z2Z3ZQ + Z1Z4ZG + Z2ZiZQ + Z3Z4ZQ + Z1Z5ZG + Z2Z5ZQ + Z3Z5ZQ (139) 
Finally using the saddle point 

1 



' 6 

we find 

C ' 



(140) 

(2) ^ 

A related, but technically different problem arise in computing series of the form 

E i-^f^'^N, (142) 

Af=l,oo 

when D is large and Aj^ are defined in eq. 11351 As soon as r > 1 the series is divergent for large 
A^. On the other end the alternating sign structure implies a certain amount of cancellation. 
This kind of computation can be found in [12], we report here for completeness. 

At this end it is useful to observe that using the convolution theorem, starting from the 
definition in that eq. (1135P we find that 

An = j d^'plipf, (143) 

where I{p) is the Fourier transform of the measure dv{x). We thus find 

G = j d°p{l + r^I{p))-^ . (144) 

After a short computation we find that I{p) is proportional to a Bessel function. It is positive 
at small p and it has a negative minimum at p = Pm, where Pm behaves as 

Pm = AD + Oil), A=^ (145) 

The value of the function I{pm) can be estimated to be 

I{p^) = -B'^, B=y^« 1.166 (146) 

This result for the function I{p) can be obtained directly from the known formulae on Bessel 
functions. It may be instructive to derive it from scratch. 

We can start from the representation for D-dimensional Fourier transform in polar coordi- 
nates 

I{p)=j dx(l-x2)(-^"^)/%xp(zpx). (147) 
Alternatively using the same techniques as before we can write 



In both cases one finds a purely imaginary saddle point in the variable x (or A) that gives as 
usual a purely real contribution. For example the value of x at the saddle point is given by 

X = (149) 

1 + VI - 4^2 ^ ' 
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where 

Dz=p (150) 
Therefore for positive p < DA we finds that 

I{p) = eM-Dgiz)), (151) 

where the function g is given by 



1 lWr^4i2 2z^ 
= 2 2 ) - l + Vl-4.^ - ^'''^ 

On the contrary as soon as p > DA the saddle point is no more purely imaginary and the result 
is complex. In the region where p = DA + t, t being of order one when D ^ oo one finds 

lip) = eM-DgiA)hit)), (153) 

where h{t) is an oscillating function of order 1. 

The chain approximation does not present difficulties for r < B, but it diverges as soon 
r > B as the effect of the pile up of subsasyntotic terms. 
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